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I .  INTRODUCTION 


In  the  first  report  for  this  contract.  Rimer,  et  al. 
(1980)  described  the  calculations  which  carried  the  Source 
3/5  ground  coupling  calculation  for  a  one  megaton  near  sur¬ 
face  nuclear  explosion  into  the  elastic  regime.  In  this 
report  we  propagate  this  elastic  wavefield  into  the  far  field 
and  compare  these  results  with  the  ground  motion  induced  by 
only  the  airblast  part  of  the  explosion. 

Synthetic  seismograms  and  response  spectra  were  cal¬ 
culated  for  ranges  varying  from  3.5  miles  to  350  miles 
(5.63  km  to  563  km)  for  three  sites  within  Minuteman  Wings 
IV,  V  and  VI.  The  results  were  examined  for  their  dependencies 
on 


•  range 

•  crustal  structure 

•  source  (full  explosion  simulation  on  a  non¬ 
linear  soil  model  versus  the  airblast  alone 
on  an  elastic  medium) . 

The  ranges  (3.5,  10,  35,  100  and  350  miles)  and  sites 
were  suggested  by  TRW  engineers  and  approved  by  Weidlinger 
Associates.  The  earth  models  were  constructed  using  data  from 
TRW  and  from  the  literature  for  Site  3  at  Wing  V,  Site  1  at 
Wing  VI  and  Site  8  at  Wing  IV. 

Peaks  in  the  computed  pseudo  relative  velocity  (PSRV) 

spectra  fall  off  from  values  between  20  cm/sec  and  100  cm/sec 

-3  -2 

at  5.63  km  to  values  between  10  cm/sec  and  10  cm/sec  at 

563  km,  depending  on  the  earth  model.  The  rate  of  decay  with 
-2  -2  5 

range  is  about  r  to  r  .  At  5.63  km,  the  pseudo  absolute 

accelerations  (PSAA)  peak  at  0.5  g  to  5.C  g  and  decay  at  a  rate 
-3  -5-4 

of  about  r  to  10  g  to  10  g  at  563  km.  The  peaks  occur  at 
periods  near  1  second  for  the  closest  range,  near  3  seconds 
at  intermediate  distances  and  near  10  seconds  at  the  farthest 
position. 


Except  at  the  closest  position,  the  peak  amplitudes  in 
PSRV  and  PSAA  for  the  full  simulation  were  within  a  factor  of 
two  of  those  for  the  airblast  alone.  The  periods  of  maximum 
PSRV  and  PSAA  agreed  closely  for  both  sources.  The  difference 
in  close  are  at  short  periods  and  are  most  likely  due  to  the 
difference  in  response  of  the  elastic  solid  for  the  airblast 
alone  and  the  nonlinear  soil  response  in  the  full  simulation.  As 
did  Murphy  and  Bennett  (1979),  we  conclude  that  the  airblast 
alone  on  an  elastic  halfspace  provides  an  inexpensive  yet 
accurate  way  of  estimating  far  field  ground  motions  from  a 
surface  burst. 

The  report  is  divided  into  three  remaining  sections. 
Section  II  describes  the  mathematical  and  numerical  methods 
used  for  calculating  far  field  ground  motions.  The  earth 
models  (elastic  and  anelastic  crustal  models)  are  discussed 
in  the  third  section.  The  synthetic  seismograms  and  response 
spectra  are  shown  and  compared  in  the  final  section. 
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II.  COMPUTATIONAL  METHODS 


The  purpose  of  this  work  is  to  estimate  ground  motions 
from  an  explosion  on  the  earth's  surface.  To  do  this  requires 
an  earth  model  and  a  method  for  computing  displacements  for 
that  model.  This  necessarily  involves  assumptions  that  make 
the  calculations  tractable  without  sacrificing  features  that 
affect  the  observed  displacements.  Our  primary  assumptions 
are : 

1.  that  the  earth  can  be  modeled  as  a  plane¬ 
layered  half space, 

2.  that  the  dominant  ground  motions  are  surface 
wave  motions ;  and 

3.  that  we  can  restrict  our  attention  to  periods 
longer  than  one  second. 

Within  a  geologic  province,  the  plane-layered  earth 
assumption  is  fairly  good  at  long  periods.  For  short  periods, 
however,  dipping  and  curving  interfaces  and  other  irregularitie 
are  important.  Since  we  have  no  accurate  knowledge  of  these 
deviations  from  plane  layers  (or  inexpensive  means  for  in¬ 
cluding  their  effects) ,  it  is  unrealistic  to  compute  motions 
for  short  periods.  We  have  chosen  one  second  as  a  cutoff. 

Numerous  observations  and  calculations  have  shown  that 
for  these  periods  the  dominant  ground  motions  at  more  than  a 
few  kilometers  from  a  surface  source  are  surface  waves.  That 
is  not  to  say  that  the  body  waves  would  not  be  evident,  but 
that  they  may  be  safely  ignored  in  estimating  peak  ground 
motions.  Our  calculations  then  are  for  surface  waves  only, 
but  will  include  higher  modes.  The  number  of  modes  (funda¬ 
mental  plus  five  overtones)  was  chosen  to  represent  the  entire 
surface  wave  solution  for  periods  greater  than  one  second,  in 
the  next  section,  response  spectra  will  be  shown  for  shorter 
periods  because  the  calculation  for  individual  modes  was  com¬ 
puted  to  frequencies  high  enough  to  prevent  aliasing. 


Surface  waves  are  the  part  of  the  solution  obtained  by 
evaluation  of  residues.  That  is,  when  evaluating  the  wave 
number  integrals  obtained  by  spatial  and  temporal  Fourier  trans¬ 
forms,  one  finds  only  the  contribution  due  to  poles  in  the 
complex  plane.  This  is  not  the  complete  solution  as  branch 
line  integrals  contribute  also.  The  branch  line  integrals  are 
responsible  for  the  body  waves  but  also  cancel  an  acausal  long 
period  contribution  from  the  residues.  Swanger  and  Boore  (1978) 
and  Herrmann  (1977)  show  that  at  smaller  ranges  the  acausal  sur¬ 
face  wave  signal  can  be  seen  but  is  small  and  has  little  effect 
on  peak  motions.  At  large  ranges,  the  acausal  signal  is  insigni¬ 
ficant.  Because  of  this  acausal  residue  solution  in  our  surface 
wave  calculations,  the  time  series  (seismograms)  are  shifted  by 
40  seconds  in  the  displays  in  Section  IV. 

In  the  remainder  of  this  section  we  describe  the  methods 
for  computing  multimode  surface  waves  from  both  the  Source  3/5 
simulation  and  the  airblast  alone. 

2.1  PROPAGATION  OF  SURFACE  WAVES  FROM  NUMERICAL  SIMULATIONS 

The  details  for  propagating  the  elastic  wavefield  from 
a  surface  in  a  numerical  grid  is  described  in  detail  in  Bache, 
et  al.  (1980).  An  excerpt  from  their  work  is  included  in 
Appendix  A.  The  method  involves  spatial  and  temporal  convolu¬ 
tions  of  monitored  displacements  and  stresses  with  a  particular 
set  of  Green's  functions.  These  Green's  functions  are  the  dis¬ 
placements  and  stresses  due  to  forces  applied  at  the  point  to 
which  the  wavefield  is  to  be  propagated.  We  compute  only 
the  surface  wave  contribution  to  the  Green's  functions  since 
we  have  chosen  to  use  only  these  types  of  waves  in  our 
analysis . 

The  explosion  simulation  in  Rimer,  et  al.  (1980)  was 
done  in  a  material  with  constant  elastic  properties.  The 
surface  on  which  the  wave  motions  were  monitored  extended  to 


* 


f 


a  depth  of  578  meters  (Figure  1) .  Since  the  crustal  models 
for  the  sites  of  interest  do  not  have  a  single  constant  layer 
in  the  upper  578  meters,  a  means  for  making  the  transition 
from  the  near-source  region  to  the  remainder  of  the  propaga¬ 
tion  path  is  needed.  Alewine  (1974)  gives  a  transmission 
coefficient  T(ui)  across  this  boundary 


T(oj) 


h  \) 

\Cr  \l 


1/2 


/ 


where  Cg  and  CR  are  the  phase  velocities  in  the  source  and 
receiver  parts  of  the  path,  and  ft-  and  A.,  are  the  source- 

KC 

depth  independent  part  of  the  surface  wave  excitation  (Figure 
2) .  Following  Bache,  Rodi  and  Harkrider  (1978) ,  the  solution 
for  the  mixed  path,  for  example  Gr  from  Equation  (A. 11),  is 


x(z)  Jjlkr)  h^2) 


/“* s  .  “M/(i,rrs  +  >sV 

V7 


For  calculations  involving  the  mixed  path,  we  used  a  source 
region  Rg  which  was  25  percent  of  the  total  path  at  the  two 
close  ranges  (5.63  and  16.1  km)  and  was  25  km  for  the  farther 
ranges  (56.3,  161  and  563  km). 


2.2  PROPAGATION  OF  SURFACE  WAVES  FROM  AN  AIRBLAST  ON  AN 

ELASTIC  HALFSPACE 

An  aspect  of  this  study  is  to  assess  the  degree  to  which 
the  ground  motions  from  the  full  nonlinear  simulation  of  a 
large  nuclear  explosion  can  be  approximated  by  those  from  the 
airblast  alone  acting  on  a  linear  elastic  half space.  In  this 
subsection  we  show  how  surface  waves  due  to  the  airblast  over¬ 
pressure  of  Needham,  et  al.  (1975)  (used  in  the  Source  3/5  cal¬ 
culation)  are  computed. 
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Figure  2.  A  schematic  of  the  mixed  propagation  path  is  shown 


We  use  the  results  of  Bache,  Masso  and  Mason  (1977) . 
Denote  the  alrblast  overpressure  by  P(r,t),  and  define 


P< 


w)  ■  f  dt  e~iut  f  P  (r, 

*  r\  *  A 


t)  Jq (kr)  rdr 


The  vertical  surface  wave  displacement  is  then 


W  -  -in  Ar  hJ2)  (kR)  P(w) 


and  the  radial  displacement  is 


u 


e0 


h|2)  (kR) 
H(q2)  (kR) 


Anelastic  attenuation  and  mixed  path  effects  are  included  as 
in  Section  2.1. 
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III.  EARTH  MODELS 


In  this  section  we  describe  the  specific  crustal  models 
used  in  our  calculations.  The  models  are  appropriate  to  three 
locations: 


1.  Wing  V,  Site  3  (referred  to  as  W5S3) 

2.  Wing  IV,  Site  8  (W4S8) ,  and 

3.  Wing  VI,  Site  1  (W6S1) . 

The  upper  crust  at  the  three  sites  were  provided  by  TRW 
engineers.  The  lower  crust  beneath  Wing  V  is  the  Northern 
Colorado  Plateau  model  of  Keller,  et  a_l.  (1976)  .  The  North 
Dakota  crust  at  Wing  VI  was  based  on  refraction  surveys  in 
southwestern  Manitoba  and  southeastern  Saskatchewan  by  Green, 
et  al.  (1980) .  The  refraction  studies  by  Healey  and  Warren 
(1969)  in  northwestern  Missouri  provided  a  model  of  the  crust 
at  Wing  IV.  The  shear  wave  dependence  on  depth  for  the  crust 
is  illustrated  in  Figure  3  and  in  Figure  4,  where  the  upper 
10  km  is  emphasized.  All  elastic  properties  are  listed  in 
Table  1. 

Geophysical  studies  generally  yield  good  estimates  for 
elastic  properties  of  the  crust.  Studies  to  determine  the 
anelastic  attenuation  factors  y  (and  ultimately  Q)  have  met 
with  mixed  success.  At  long  periods  the  results  seem  fairly 
good.  At  periods  less  than  10  seconds,  however,  the  results 
are  scattered  with  a  corresponding  lack  of  confidence  in  their 
validity.  The  primary  difficulty  is  that  in  order  to  isolate 
the  effects  of  anelastic  attenuation  it  is  necessary  to  have 
a  thorough  understanding  of  the  elastic  properties  and  an 
ability  to  remove  their  effects.  The  elastic  effects  are 
usually  accounted  for  by  using  plane-layered  models,  so  that 
y  and  Q  end  up  accounting  not  only  for  anelastic  absorption 
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Figure  4.  The  shear  wave  dependence  on  depth  in  the 
upper  10  km  is  shown  for  the  three  sites. 
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DEPTH  (KM) 


TABLE  1 


EARTH  MODELS 

Layer 

Thickness  a  8  P  q 

(km)  (km/sec)  (km/sec)  (gm/cm3)  w6 


W5S3 


0.13 

1.10 

0.64 

1.80 

20 

0.63 

2.30 

0.95 

2.30 

20 

1.04 

3.20 

1.50 

2.60 

35 

0.50 

4.40 

2.30 

2.65 

50 

2.70 

5.50 

2.90 

2.70 

250 

4.00 

6.10 

3.50 

2.75 

400 

15.00 

6.40 

3.70 

2.80 

1200 

16.00 

6.80 

3.90 

2.90 

2000 

1.00 

7.80 

4.50 

3.20 

2000 

W6S1 


0.24 

1.70 

0.74 

2.40 

20 

0.22 

2.20 

1.19 

2.60 

20 

0.30 

3.20 

1.70 

2.70 

35 

0.24 

5.50 

3.20 

2.70 

225 

1.00 

6.00 

3.50 

2.70 

300 

7.00 

6.20 

3.60 

2.80 

400 

16.00 

6.60 

3.80 

2.90 

1200 

16.00 

6.60 

3.80 

2.90 

1800 

1.00 

8.00 

4.60 

3.20 

2000 

W4S8 


0.053 

3.  51 

1.77 

2.50 

20 

0.221 

3.93 

2.10 

2.67 

20 

0.  086 

4.  25 

2.27 

2.67 

20 

0.240 

4.68 

2.50 

2.67 

20 

2.  000 

5.60 

3.23 

2.68 

50 

2.  400 

5.60 

3.23 

2.68 

250 

5.  000 

6.10 

3.50 

2.69 

300 

10. 000 

6.20 

3.  60 

2.70 

1200 

20.000 

6.60 

3.81 

2.90 

1800 

1.000 

8.00 

4.62 

3.30 

2000 

but  also  for  elastic  scattering  and  curving  or  dipping  inter¬ 
faces.  Because  our  choice  of  Q  model  strongly  affects  seismic 
amplitudes  it  is  important  that  we  make  the  best  choice  possi¬ 
ble. 

Although  this  is  still  a  relatively  unexplored  field, 
especially  at  short  periods,  we  can  use  the  results  of  Bache, 
Swanger  and  Shkoller  (1980) .  In  their  work,  they  compared 
regional  synthetic  and  observed  surface  wave  recordings  for 
periods  around  one  second.  They  deduce  a  Q  model  based  on 
matching  of  frequency  content,  peak  amplitude,  and  signal 
duration.  The  Q  models  used  in  our  calculations  are  based 
on  techniques  developed  in  their  work.  Figure  5  shows  the  Q 
models. 

A  shortcoming  of  the  synthetic  seismograms  computed  by 
Bache,  Swanger  and  Shkoller  (1980)  is  that  the  signal  duration 
is  shorter  than  the  observations.  Presumably  the  seismograms 
we  have  computed  for  this  project  suffer  in  the  same  way. 
However,  we  can  show  that  this  has  little  affect  on  the  re¬ 
sulting  response  spectra.  Consider  the  seismogram  shown  in 
Figure  6  which  is  a  recording  of  the  SALMON  explosion  at  a 
range  of  244  km.  This  recording  was  made  by  a  short  period 
instrument  whose  response  peaks  at  3  Hz  which  greatly  attenu¬ 
ates  the  longer  period  dominant  surface  waves  and  enhances  the 
body  wavos.  The  surface  wave  portion  of  this  record  begins 
roughly  with  signals  arriving  at  an  apparent  velocity  of  3.5 
km/sec  (70  seconds) .  Synthetic  surface  wave  seismograms  show 
ground  motion  dominated  by  waves  traveling  with  speeds  in  the 
range  3.5  to  2.9  km/sec  (84  seconds).  The  observed  ground 
motions  last  much  longer.  Response  spectra  of  the  time  inter¬ 
val  from  70  to  84  seconds  and  for  a  section  twice  as  long  from 
70  to  98  seconds  are  shown  in  Figures  7  and  8.  Clearly  the 
two  spectra  are  very  close  except  at  the  periods  longer  than 
the  sample  intervals.  Thus  the  weakness  that  calculated 
seismograms  using  the  Q  model  of  Bache,  Swanger  and  Shkoller 


■*  -Wv  « 


DEPTH  (KM) 


W553 

W6S1 

W4S8 


Q 


0.0  400.0  800.0  1200.0  1600.0  2000  0 


Q 


Figure  5.  The  anelastic  attenuation  parameter  Q  is 
plotted  versus  depth. 
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Figure  8.  The  response  spectrum  of  the  portion  of 
the  seismogram  in  Figure  6  between  70 
and  98  seconds  is  shown. 


(1980)  do  not  match  the  observed  signal  durations  does  not 
affect  the  peak  period  and  amplitude  of  the  corresponding 
response  spectra. 

The  relative  values  of  Q  in  these  models  are  more 
closely  constrained  than  the  absolute  values.  To  estimate 
the  effect  of  having  chosen  our  values  of  Q  to  be  a  factor 
of  two  too  large  (an  extreme) ,  we  have  computed  seismograms 
using  our  Q  model  scaled  by  one-half. 


IV.  GROUND  MOTION  CALCULATIONS 


The  objective  of  the  calculations  described  in  this 
section  is  to  determine  the  differences  in  ground  motion  due 
to 


•  the  full  simulation  of  the  explosion  versus 
the  airblast  on  an  elastic  medium, 

•  the  distance  from  ground  zero,  and 

•  the  crustal  models  that  typify  the  Minuteman 
wings  chosen  for  this  study. 

We  begin  our  discussion  with  an  analysis  of  the  ground  motion 
from  the  full  surface  burst  simulation  versus  that  produced 
by  only  the  airblast. 

4.1  GROUND  MOTIONS  FROM  THE  FULL  SIMULATION  AND  THE 

AIRBLAST  OVER  AN  ELASTIC  MEDIUM 

This  comparison  is  performed  on  the  W5S3  structure. 
Figures  9  and  10  show  the  calculated  vertical  displacements 
at  the  five  ranges  for  the  two  sources.  Figures  11a  through 
lie  compare  the  corresponding  response  spectra.  The  airblast 
alone  calculation  was  done  with  exactly  the  sarnie  mixed  propa¬ 
gation  path  and  overpressure  function  (from  Needham,  et  al . , 
1976)  that  were  used  for  the  full  simulation.  The  differences 
in  response  are  due  only  to  the  source  (airblast  alone  versus 
airblast  plus  radiation  effects)  and  to  the  soil  response 
(elastic  versus  nonlinear  near  source  soil  models) . 

The  peak  displacements  from  the  seismograms  for  the  full 
simulation  range  from  a  few  centimeters  at  r  *  5.63  km  to  a 
few  microns  at  r  **  563  km.  Those  from  the  airblast/elastic 
source  are  the  same  at  r  =  56.3  and  two  to  four  times  as  great 
at  the  other  ranges.  The  shapes  of  the  seismograms  are  basi¬ 
cally  very  similar  but  have  many  differences  in  detail. 

Murphy  and  Bennett  (1979)  made  similar  comparisons 
between  the  vertical  displacements  at  1800  meters  from  another 
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Vertical  displacement  (microns)  are  shown  for  the  full  Simula 
tion  at  W5S3  for  five  ranges.  Tick  marks  are  at  40  second 
intervals,  beginning  at  -40  seconds. 
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Figure  10.  Vertical  displacements  (microns)  for  the  airblast  alone  on  an 
elastic  medium  with  mixed  propagation  path  are  shown  for  five 
ranges.  Tick  marks  are  at  40  second  intervals,  beginning  at 
-40  seconds. 


Figure  11a.  Response  spectra  for  the  full  simulation  (left)  is  compared 
with  the  airblast  alone  solutions  for  the  W5S3  earth  model 


Figure  lid.  Response  spectra  for  the  full  simulation  (left)  is  compared 
with  the  airblast  alone  solutions  for  the  W5S3  earth  model 
at  r  =  161.0  km. 


Figure  lie.  Response  spectra  for  the  full  simulation  (left)  is  compared  with 
the  airblast  alone  solutions  for  the  W5S3  earth  model  at  r  * 


nonlinear  finite  difference  simulation  of  the  airblast  loading 
(by  Sandler,  1978)  and  from  an  airblast  over  an  elastic  half¬ 
space.  They  found  that  the  general  character  of  the  ground 
motions  from  both  sources  was  in  agreement  and  that  the  peak 
displacements  for  the  elastic  case  were  two  to  three  times 
those  for  the  nonlinear  case. 

The  response  spectra  is  compared  in  Figures  11a  through 
lie.  Their  shapes  are  quite  similar  at  the  three  farthest  ranges 
and  at  long  periods  for  the  short  ranges.  The  short  period 
response  at  ranges  5.63  and  16.1  km  are  very  different.  One 
expects  that  the  differences  between  the  linear  and  nonlinear 
soil  responses  to  be  greatest  at  high  frequencies.  Due  to  the 
anelastic  absorption  operator  e~Yr,  the  higher  frequency 
energy  is  attenuated  more  with  range  which  therefore  suppresses 
the  high  frequency  differences. 

The  peak  amplitudes  and  corresponding  periods  in  the 
response  spectra  are  tabulated  in  Table  2  and  plotted  in 
Figures  12a  through  12c.  Generally,  the  peaks  for  the  two 
source  models  occur  at  the  same  periods.  The  amplitudes  are 
generally  one  and  one-half  to  two  times  greater  for  the  air- 
blast/elastic  case  than  for  the  full  simulation.  Thus  the 
nonlinear  absorption  at  the  soils  beneath  the  explosion  cause 
the  coupling  into  seismic  energy  to  be  reduced  by  about  a 
factor  of  two  without  shifting  the  periods  at  which  peak 
motions  occur,  except  at  the  closest  distance. 


4.2  THE  MIXED  PROPAGATION  PATH 

Recall  that  in  Section  II  we  introduced  the  mixed 
propagation  path  to  obtain  a  transition  from  the  source 
region  to  the  remainder  of  the  propagation  path.  The  mix  of 
source  and  receiver  path  was  chosen  arbitrarily.  We  can 
check  the  effects  of  our  choice  by  comparing  the  motions  from 
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TABLE  2 


PERIODS  AND  AMPLITUDES  AT  PEAKS  IN  VERTICAL  RESPONSE 
SPECTRUM  FOR  THE  FULL  SIMULATION  AND  THE 
AIRBLAST/ELASTIC  CASE  FOR  THE  W5S3  EARTH  MODEL 


FULL  SIMULATION  AIRBLAST/ELASTIC 

PRD 


Range 

T 

max 

Amax 

(cm) 

T 

max 

Amax 

(cm) 

5.63 

2.5 

6.0 

2.2 

11.2 

16.1 

2.8 

1.86 

2.8 

2.7 

56.3 

4.0 

1.4 

x  10-1 

4.0 

1.6 

x  10"1 

161.0 

5.0 

1.6 

-2 

x  10  z 

5.0 

2.2 

x  10"2 

563.0 

5.6 

1.3 

x  10"3 

6.3 

2.5 

x  10”3 

PSRV 

Range 

T 

max 

A 

max 

(cm/ sec) 

T 

max 

\iax 

(cm/sec) 

5.63 

2.0 

18.0 

2.0 

34.3 

16.1 

2.8 

4.1 

2.2 

6.8 

56.3 

4.0 

2.2 

x  10"1 

4.0 

2.5 

x  10"1 

161.0 

5.0 

2.0 

x  10"2 

5.0 

2.8 

x  10-2 

563.0 

5.6 

1.5 

x  10"3 

5.6 

2.8 

x  10-3 

PSAA 

Range 

T 

max 

A 

max 

<g> 

T 

max 

Amax 

(g> 

5.63 

2.0 

5.8 

x  10"2 

1.3 

4.9 

x  10-2 

16.1 

2.2 

1.1 

x  10"2 

2.0 

1.9 

x  10-2 

56.3 

4.0 

3.5 

x  10-4 

4.0 

4.1 

x  10"4 

161.0 

5.0 

2.6 

x  10"5 

5.0 

3.6 

x  10-5 

563.0 

5.6 

1.7 

x  10'6 

5.6 

3.2 

x  10’6 

Figure  12a.  Peak  values  in  PRD  are  plotted  versus 
range  for  the  full  simulation  and  the 
airblast/elastic  case.  In  the  lower 
plot,  the  values  are  normalized  to  those 
of  the  full  simulation. 
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Figure  12b.  Peak  values  of  PSRV  are  plotted  versus 
range  for  the  full  simulation  and  the 
airblast/elastic  case.  In  the  lower 
plot,  the  values  are  normalized  to  those 
of  the  full  simulation. 
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the  airblast/elastic  source  with  and  without  the  mixed  path. 
That  is,  we  compare  the  above  results  discussed  with  those 
for  which  the  source  region  radius  Rg  goes  to  zero.  The 
medium  is  then  a  simple  layered  halfspace  and  no  transition 
is  needed. 

We  again  use  the  W5S3  structure.  Figure  13  shows  the 
synthetic  seismograms  for  the  single  path.  These  should  be 
compared  with  the  seismograms  in  Figure  10.  The  response 
spectra  are  compared  in  Figures  14a  through  14e  and  peaks  are 
plotted  versus  range  in  Figures  15a  through  15c.  The  seismo¬ 
gram  shapes  are  similar  and,  except  at  the  smallest  distance, 
the  peak  displacements  are  close  to  those  obtained  from  the 
seismograms  in  Figure  10.  The  peak  response  spectra  agree 
well  except  for  the  PSAA  at  5.63  km.  This  peak  occurs  below 
the  one  second  cutoff  and  is  not  reliable.  If,  instead,  we 
use  the  value  at  1.2  seconds  where  the  mixed  path  spectrum 
peaked,  the  ratio  of  the  single  path  to  mixed  path  peaks  is 
about  three.  Since  the  path  is  mixed  only  in  the  top  600  m, 
we  expect  its  effects  to  be  seen  only  at  high  frequencies  and 
small  ranges  which  is  the  case.  So  except  for  the  accelera¬ 
tions  at  r  =  5.63,  our  choice  of  source  region  dimension  is 
not  important  and  for  the  purpose  of  computing  ground  motions 
for  the  airblast/elastic  case  we  can  use  the  simpler  single 
path  model. 

4.3  COMPARISON  OF  GROUND  MOTIONS  FOR  THREE  EARTH  STRUCTURES 

In  this  section  we  explore  the  differences  in  response 
of  the  three  crustal  structures  discussed  in  Section  II 
(Figures  3  and  4) .  For  these  calculations  we  have  used  the 
airblast/elastic  source  and  the  single  propagation  path. 

Note  that  the  structures  are  very  close  to  being  the 
same  below  about  10  km  (Figure  3)  so  that  differences  in 
response  must  be  due  largely  to  differences  above  that  depth. 
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Figure  14a.  Response  spectra  for  the  airblast  on  an  elastic  medium  for 
the  mixed  path  (left)  and  for  the  single  path  are  shown; 
r  =  5.63  km. 
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Figure  15a.  Peak  values  of  PRD  are  plotted  versus 
range  for  the  airblast/elastic  source 
with  a  single  and  a  mixed  propagation 
path.  Values  in  bottom  plot  are 
normalized  to  the  mixed  path  values. 
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Figure  15b.  Peak  values  of  PSRV  are  plotted  versus 
range  for  the  airblast/elastic  source 
with  a  single  and  a  mixed  propagation 
path.  Values  in  bottom  plot  are 
normalized  to  the  mixed  path  values. 
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Figure  15c.  Peak  values  of  PSAA  are  plotted  versus 
range  for  the  airblast/elastic  source 
with  a  single  and  a  mixed  propagation 
path.  Values  in  bottom  plot  are 
normalized  to  the  mixed  path  values. 


At  W5S3  and  W6S1/  the  shear  wave  speeds  (6)  are  very  low  at 
the  surface  (0.6  to  0.7  km/sec).  W5S3  has  a  relatively  gradual 
increase  in  8  compared  to  W6S1  which  quadruples  its  value  with¬ 
in  1  km  of  the  surface.  At  W4S8,  the  surface  value  of  8  is 
1.8  km/sec,  increases  to  3.2  km/sec  in  the  upper  1  km  and 
gradually  increases  below  that. 

The  most  obvious  differences  in  the  synthetic  seismo¬ 
grams  (Figures  16  through  18)  is  the  signal  duration  at  W4S8 
compared  to  W6S1  and  W5S3.  At  all  ranges,  the  signal  at  W4S8 
is  rather  pulse  like  without  the  later  time  ringing  evident 
at  the  other  sites.  The  ringing  is  due  to  waves  traveling 
in  the  slow  upper  1  km  of  crust  at  W5S3  and  W6S1. 

Differences  in  excitation  of  the  crustal  structure  are 
very  apparent  in  the  response  spectra  compared  in  Figures  19a 
through  19e.  Amplitudes  and  periods  of  spectral  peaks  are 
given  in  Table  3  and  amplitudes  are  plotted  in  Figures  20a 
through  20c.  For  the  most  part,  the  response  spectra  peak 
in  the  period  range  of  interest  (one  to  ten  seconds) .  Peaks 
in  the  response  spectra  at  W5S3  occur  at  longer  periods  than 
at  the  other  wings,  while  the  peaks  at  W6S1  and  W4S8  occur  at 
roughly  the  same  period  (except  at  the  closest  receiver) .  The 
dominance  of  energy  at  longer  periods  at  VJ5S3  is  due  to  the 
gradual  increase  in  8  in  the  upper  6  km  of  the  crust  relative 
to  the  other  sites. 

From  Figures  20a  through  20c,  we  see  a  tendency  for 
peak  spectral  amplitudes  at  W4S8  and  W6S1  to  increase  relative 
to  W5S3  as  a  function  of  range.  The  values  of  peak  PSRV  at 
W4S8  and  W6S1  are  within  a  factor  of  two  of  those  at  W5S3. 

The  PRD  and  PSAA  peaks  differ  by  as  much  as  a  factor  of  four, 
with  the  greatest  difference  in  close  for  the  PRD  peaks  and  at 
the  greater  ranges  for  the  PSAA. 


Figure  16.  Synthetic  displacement  seismograms  at  W5S3  for  the  airblast/ 
elastic  case  with  a  single  propagation  path  are  shown.  Tick 
marks  are  at  40  second  intervals,  beginning  at  -40  seconds. 


Figure  17.  Synthetic  displacement  seismograms  at  W4S8  for  the  airblast/ 
elastic  case  with  a  single  propagation  path  are  shown.  Tick 
marks  are  at  40  second  intervals,  beginning  at  -40  seconds. 
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TABLE  3 


PERIODS  AND  AMPLITUDES  AT  PEAKS  IN  VERTICAL  RESPONSE 
SPECTRA  FOR  THE  AIRBLAST-ALONE  SOURCE  AT  THE 
THREE  WINGS  USING  SINGLE  PATH  PROPAGATION 


PRD 

W5S3 

W4S8 

W6S1 

Range 

T  A 

max  max 

Tmax  ^max 

T  A 

max  max 

5.63 

0.9 

1.0 

X 

105 

0.5 

2.6 

X 

104 

1.0 

8.3 

X 

104 

16.1 

2.5 

1.8 

X 

104 

0.7 

5.2 

X 

103 

1.3 

1.1 

X 

104 

56.3 

2.8 

1.4 

X 

103 

2.0 

7.7 

X 

IQ2 

1.7 

7.8 

X 

102 

161.0 

5.6 

2.1 

X 

102 

3.6 

2.2 

X 

102 

5.6 

2.0 

X 

102 

563.0 

10.0 

4.1 

X 

101 

7.1 

7.3 

X 

IO1 

6.3 

7.6 

X 

IO1 

PSRV 

5.63 

0.9 

7.2 

X 

105 

0.5 

3.7 

X 

105 

1.0 

5.2 

X 

105 

16.1 

2.2 

4.6 

X 

io4 

0.7 

4.6 

X 

104 

1.3 

5.4 

X 

104 

56.3 

2.8 

3.0 

X 

103 

1.1 

3.3 

X 

103 

1.6 

3.1 

X 

103 

161.0 

4.5 

2.5 

X 

102 

2.5 

4.8 

X 

102 

3.2 

2.7 

X 

102 

563.0 

8.9 

2.8 

X 

io1 

5.62 

7.0 

X 

IO1 

5.6 

8.0 

X 

IO1 

PSAA 

5.63 

0.7 

5.1 

X 

106 

0.5 

3.7 

X 

105 

0.7 

4.3 

X 

106 

16.1 

0.7 

2.5 

X 

io5 

0.5 

4.8 

X 

105 

1.3 

2.7 

X 

105 

56.3 

2.8 

6.7 

X 

io3 

1.0 

1.9 

X 

104 

1.6 

1.2 

X 

10  4 

161.0 

3.6 

3.7 

X 

102 

1.8 

1.5 

X 

IO3 

3.2 

5.3 

X 

102 

563.0 

7.1 

2.3 

X 

101 

4.5 

8.5 

X 

IO1 

4.5 

9.7 

X 

IO1 
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Figure  19a.  Response  spectra  for  the  airblast/elastic  case  with 
propagation  path  are  shown  for  the  three  sites  at  r 


Figure  19c.  Response  spectra  for  the  airblast/elastic  case  with 
propagation  path  are  shown  for  the  three  sites  at  r 


Figure  19d.  Response  spectra  for  the  airblast/elastic  case  with 
propagation  path  are  shown  for  the  three  sites  at  r 


6S1 


Figure  19e.  Response  spectra  for  the  airblast/elastic  case  with  a  sing 
propagation  path  are  shown  for  the  three  sites  at  r  =  563. 
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Figure  20a.  Peak  values  of  PRD  are  plotted  versus 
range  for  W5S3,  W4S8  and  W6S1.  Values 
are  normalized  to  W5S3  in  the  bottom 
plot. 
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Figure  20b. 


Peak  values  of  PSRV  are  plotted  versus 
range  for  W5S3,  W4S8  and  W6S1.  Values 
are  normalized  to  W5S3  in  the  bottom 
ot. 
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Figure  20c. 


Peak  values  of  PSAA  are  plotted  versus 
range  for  W5S3,  W4S8  and  W6S1.  Values 
are  normalized  to  W5S3  in  the  bottom 
plot. 
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4.4  DEPENDENCE  ON  Q 

The  Q  models  described  in  Section  III  were  derived  from 
a  body  of  evidence  which  constrains  the  relative  values  of  Q 
within  the  stack  of  layers  considerably  more  than  their  abso¬ 
lute  values.  To  provide  a  bound  on  the  greatest  likely  attenua¬ 
tion  we  have  computed  seismograms  for  the  W5S3  crustal  structure 
in  which  all  values  of  Q  have  been  halved  (or  equivalently,  the 
values  of  y  doubled) .  The  airblast/elastic  source  and  the 
single  path  propagation  model  have  been  used. 

A  list  comparing  amplitudes  and  periods  of  peaks  in  the 
response  spectra  for  the  Q  model  in  Table  1  and  for  the  half-Q 
model  is  given  in  Table  4.  Peak  amplitudes  versus  range  are 
shown  in  Figures  21a  through  21c. 

The  periods  T _ _  shift  very  little  for  the  PRD  and  PSRV 

max 

for  all  ranges  and  for  the  PSAA  at  the  closer  ranges.  The 

peaks  are  attenuated  at  the  two  farthest  stations  by  factors 

ranging  from  1.3  to  2.6,  and  by  factors  as  large  as  5.6  at 

the  three  inner  stations.  The  inner  ranges  are  affected  more 

because  the  parameter  y  is  a  rapidly  increasing  function  of 

frequency.  For  example,  the  value  of  y  for  the  fundamental 

-4 

mode  at  8.9  seconds  is  5.1  x  10  km  while  at  0.9  seconds,  y 
is  1.7  x  10  1  km.  Thus  at  r  =  5.63,  where  the  response 
spectra  peaks  at  0.9  seconds,  the  attenuation  due  to  doubling 
Y  is  about  e-yr /e~^r ,  or  0.38.  At  r  *  5.63  km  on  the  other 
hand  the  peak  is  at  8.9  seconds  and  is  reduced  by  0.75,  about 
half  the  reduction  seen  at  5.63  km. 

It  is  well  known  that  the  Q  model  is  a  critical  param¬ 
eter  required  for  an  accurate  estimate  of  seismic  ground 
motion.  Recently  new  techniques  have  become  available  for 
determining  Q  on  a  site  specific  basis.  Consideration  should 
be  given  to  their  application  to  those  sites  where  ground 
motion  predictions  of  the  type  discussed  here  are  required. 
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TABLE  4 


PERIODS  AND  AMPLITUDES  AT  PEAKS  IN  VERTICAL  RESPONSE 
SPECTRA  FOR  THE  AIRBLAST  ALONE  SOURCE  AT  WING  5 
FOR  THE  Q  MODEL  IN  FIGURE  5  AND  FOR 
THAT  MODEL  WITH  ALL  VALUES  OF  Q  DIVIDED  BY  TWO 


Range 


T 

max 


Q 

A 

max 


PRD 

HALF-Q 

T  A 

max  max 


Amax  (0)"WHal£-Q> 


5.63 

0.9 

1.0 

X 

105 

0.9 

3.0 

X 

104 

3.3 

16.1 

2.5 

1.8 

X 

104 

2.5 

7.2 

X 

103 

2.5 

56.3 

2.8 

1.4 

X 

103 

3.2 

3.5 

X 

102 

4.0 

161.0 

5.6 

2.1 

X 

102 

5.6 

1.1 

X 

102 

1.9 

563.0 

10.0 

4.1 

X 

101 

10.0 

3.1 

X 

101 

1.3 

PSRV 

5.63 

0.9 

7.2 

X 

105 

0.9 

2.1 

X 

5 

10 

3.4 

16.1 

2.2 

4.6 

X 

104 

2.5 

1.8 

X 

104 

3.3 

56.3 

2.8 

3.0 

X 

103 

3.2 

7.0 

X 

102 

4.3 

161.0 

4.5 

2.5 

X 

102 

5.6 

1.3 

X 

102 

1.9 

563.0 

8.9 

2.8 

X 

101 

8.9 

1.9 

X 

101 

1.5 

PSAA 

5.63 

0.7 

5.1 

X 

106 

0.9 

1.5 

X 

106 

3.4 

16.1 

0.7 

2.5 

X 

105 

2.5 

4.5 

X 

104 

5.6 

56.3 

2.8 

6.7 

X 

103 

3.2 

1.4 

X 

103 

4.8 

161.0 

3.6 

3.7 

X 

102 

5.6 

1.4 

X 

102 

2.6 

563.0 

7.1 

2.3 

X 

101 

8.9 

1.4 

X 

101 

1.6 

56 


*  '  *^5*  **  *  -Ji  ****!*--  •***#•***'»/-;-  - 


„  LEGEND 

□  =  WING5,  AIRBLAST  ALONE.  SINGLE  PATH 
0=  WING5,  ARBL..  SNGL  PATH.  Q/2 


id  Iff 

RANGE  (KM) 


id  id 

RANGE  (KM) 


Figure  21a. 


Peak  values  of  PRD  are  plotted  for  W5S3 
using  the  Q  model  in  Figure  5  and  using 
that  model  with  all  values  of  Q  divided 
by  two. 
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Figure  21b.  Peak  values  of  PSRV  are  plotted  for  W5S3 
using  the  Q  model  in  Figure  5  and  using 
that  model  with  all  values  of  Q  divided 
by  two. 
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Figure  21c. 
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Peak  values  of  PSAA  are  plotted  for  W5S3 
using  the  Q  model  in  Figure  5  and  using 
that  model  with  all  values  of  Q  divided 
by  two . 
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APPENDIX  A 


SYNTHETIC  SEISMOGRAMS 
FROM  COMPLEX  SOURCE  CALCULATIONS 

A.l  INTRODUCTION 

There  is  increasing  interest  in  the  use  of  large 
scale  finite  difference  or  finite  element  computer  pro* 
grams  to  simulate  the  complex  nonlinear  material  behavior 
that  occurs  in  the  region  immediately  surrounding  earth¬ 
quake  or  explosion  sources.  To  fully  appreciate  the  sig¬ 
nificance  of  the  output  of  these  programs,  it  is  important 
to  be  able  to  compute  the  ground  motions  at  ranges  much 
larger  than  the  maximum  dimension  of  the  computational 
grid.  As  long  as  the  material  response  can  be  assumed  to 
be  linearly  elastic  outside  the  immediate  source  region, 
analytical  techniques  can  be  used  to  continue  the  wave 
field  to  the  ranges  of  interest.  In  this  Appendix  we 
describe  these  techniques. 

When  the  source  is  assumed  to  be  embedded  in  an 
elastic  wholespace,  an  "equivalent  elastic  source"  repre¬ 
sentation  can  be  constructed  by  expanding  the  outgoing 
wave  field  in  spherical  harmonics  (e.g.,  Bache  and  Hark- 
rider,  1976).  The  source  is  then  expressed  in  terms  of 
multipolar  coefficients  that  can  be  used  as  input  to  pro¬ 
grams  for  computing  synthetic  seismograms.  This  technique 
has  been  used  to  study  the  results  from  complex  explosion 
(Cherry,  et.  al.,  1975,  1976;  Bache  and  Masso,  1978)  and 
earthquake  (Bache,  et.  al . ,  1980)  sources  computed  with 
finite  difference  programs. 

When  boundaries  (e.g.,  a  free  surface)  are  present 
in  the  vicinity  of  the  source,  the  ground  motions  at 
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ranges  outside  the  numerical  grid  can  be  computed  with  a 
wave  field  continuation  method.  Based  on  an  elastoaynamic 
representation  theorem  (e.g.,  Burridge  and  Knopoff, 

1964) .  This  method  requires  the  convolution  of  stresses 
and  displacements  from  the  elastic  region  of  the  numerical 
source  calculation  with  Green's  functions  representing  the 
elastic  wave  propagation  through  the  medium  outside  the 
source  region. 

In  subsequent  sections  we  begin  with  the  formulation 
of  the  basic  theory.  For  earthquake  sources  the  wave 
fiela  continuation  must  be  done  in  three  dimensions.  How¬ 
ever,  many  important  complex  explosion  phenomena  are  axi- 
symmetric  and  the  theory  is  specialized  for  this  case. 

Some  of  the  development  for  the  axisymmetric  case  can  also 
be  used  for  three-dimensional  sources  in  plane-layered 
earth  models,  since  the  Green's  functions  are  azimuth-in¬ 
dependent  for  such  structures. 

The  details  of  the  implementation  depend  on  the 
choice  of  Green's  functions  for  the  propagation.  Our 
emphasis  here  is  on  the  Rayleigh  wave  modes  for  plane¬ 
layered  earth  models.  Test  problems  which  illustrate  this 
case  and  demonstrate  the  accuracy  that  can  be  achieved  are 
presented,  in  Bache,  et  al.  (1980)  . 


A.  2  FORMULATION 


For  the  analytic  continuation  we  use  an  elastody- 
namic  representation  theorem  (Burridge  and  Knopoff, 

1964).  Assuming  no  body  forces  are  present,  this  theorem 
may  be  written 


ls  [Gi*T" 
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si**ui  “ 
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( A.  1) 
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where  the  repeated  indices  imply  summation  and  the 
following  notation  has  been  used: 

uF(x  , t) :  Components  of  the  displace 

'  ivf  '  ' 

ments  at  xf  outside  the 
source  region; 


G^lx  ,x  ,t):  Displacement  in  the  direction 

3  j  at  xM  due  to  a  unit 

impulsive  force  at  Xp  in  the 
direction  i; 


Tj <*„>*>= 


n 


lc: 


The  components  of  the  stress  ten¬ 
sor  at  x„  due  to  a  unit  impulsive 

M 

force  at  Xp  in  direction  i; 

Components  of  the  monitored 
traction  vector  on  the 
monitoring  surface; 

Components  of  the  monitored 
displacement  vector  on  the 
monitoring  surface; 

Components  of  the  normal  to 
the  monitoring  surface,  with 
the  positive  sense  being  for 
vector  pointing  away  from  the 
source; 


Area  of  the  monitoring  surface. 


Also,  convolution  is  indicated  by 


f  *g 


f(t-t)g(x)  dx 


(A. 2) 
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Thus,  the  solution  at  some  distant  point  is  obtained 
by  first  convolving  suitable  Green's  functions  with 
stresses  and  displacements  monitored  on  some  surface  that 
encloses  the  region  of  inelastic  material  deformation* 

The  final  solution  is  then  obtained  by  integrating  the 
convolution  products  over  the  monitoring  surface.  The 
restrictions  on  this  technique  are,  first,  that  the  region 
outside  the  monitoring  surface  must  be  linearly  elastic. 
Second,  if  the  problem  domain  is  bounded,  the  Green's 
functions  must  satisfy  appropriate  boundary  data.  The 
elastic  medium  may  be  anisotropic  or  have  discontinuous 
properties,  though  the  computation  of  Green's  functions  is 
much  more  difficult  for  such  media  than  for  homogeneous 
isotropic  media. 

To  this  point  the  formulation  is  entirely  general. 
Let  us  now  consider  cases  where  the  source  and  propagation 
medium  are  axisymmetr ic.  This  is  the  appropriate  geometry 
for  explosions  near  the  free  surface  (or  other  horizontal 
boundaries)  in  plane-layered  earth  models. 

Stresses  and  displacements  are  monitored  on  some 
surface  with  a  typical  point  M.  We  wish  to  compute  the 
displacement  field  at  F.  Let  the  monitored  quantities  at 
M  be  denoted  by  o^_,  0**  ,  u*1  ana  u^.  Then 

the  components  of  the  traction  vector  and  displacement, 
referred  to  a  local  coordinate  system  (x^,  x^,  x^) 
at  F,  are  given  by: 

T^1  «  (o*?_  cose  +  sine)  cos  5  , 

T2  “  "  (orr  cose  +  °r 2  sin6)  sin  €  * 

*  of,  cose  +  0**  sine  , 

3  r  z  zz 
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Ur  cos  £  » 
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ut  sin  £  , 


(A. 3) 


M  M 
u3  *  uz  ' 


=  cos0  cos  £  ,n2  ■  -cose  sin  £  ,  n3  ■  sine 

Also,  the  displacements  at  F  in  the  global  R-z  system  ate 
related  to  the  displacements  in  the  x^,  x2,  x^ 
system  by 

(A. 4) 


F  F  F 

Ur  -  cos*  -  u2  sin  *  , 


Using  these  coordinate  transformations,  (A.l)  may  be  written: 

Up  *  -  f  r  |"(GR*<jM  +  GR*oM  )  -  (SR  *uM  +  SR  *uM)  1  dz 
R  J s  L  c  rr  z  rz7  v  rr  r  rz  z7  J 

M 

-  f  r  T(GR*aM  +  GR*aM  )  -  (SR  *uM  +  SR  *uM)  1  dr  , 
J K  L  *  *z  r  rz'  1  rz  r  zz  z7  J  ' 

M 

u*  ■ -/' [lG**°”  +  G**°«)  '  (S?t*u?  *  S?**V ]  dI 

-/  c  [tGI*°2Z  +  G?*°?z>  *  ‘Srz*u?  +ur  *  slz*“z>  ] 


R  R 

where  the  quantities  Gr,  Gz,  etc.,  are  azimuthal 
averages  of  the  Green's  functions  in  the  x^,  x2,  x3 
system.  The  explicit  formulas  for  the  azimuthally 
averaged  Green's  functions  are  the  following: 


dr 


X3  *  ZM  "  ZF  *  (A*8) 

For  axisymmetr ic  problems  formulated  to  compute  the  solu¬ 
tion  in  the  w ,  k  domain  (e.a.,  modal  solutions),  the  azi¬ 
muthal  integrations  in  (A. 6)  and  (A. 7)  may  be  computed 
analytically.  In  the  p,  t  domain  (e.g.,  generalized  ray 
theory),  these  integrations  cannot,  in  general,  be  per¬ 
formed  analytically.  However,  in  the  far-field,  an 
analytic  approximation  for  these  integrals  may  be  used. 

The  formulas  for  The  azimuthally  averaged  stress  Green's 
R 

functions,  S  ,  etc.,  are  given  by: 

1  2  i  1.2 

(S^  cos  5  -  sm25  +  S*2  sm  5  cos^ 

-  (S^  cos2?  -  S21  sin2?  +  S22  sin2?  )  sin<|> 
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(S^  cos?  -  Sj2  sin?)  cosi^ 

(S2^  cos?  -  Sj2  sin?)  sin'l' 


d$  ; 


(A. 9) 
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If  the  azimuthal  averages  in  (A. 6)  and  (A. 7)  can  be 
analytically  determined,  it  is  not  necessary  to  use 
(A. 9).  instead,  the  constitutive  equation  of  the  medium 
at  the  monitoring  point  can  be  used.  Thus,  for  an  iso¬ 
tropic  medium 
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(A. 10) 


68 


1' 

*\ 


where  X  and  y  are  Lame's  constants,  and  formulas  of  the 

R  R 

same  type  can  be  written  for  S  ,  and 

CR  rc  2r 

Szz  * 

The  basic  formulation  for  analytical  continuation  of 
axisymmetric  source  computations  is  now  complete.  These 
equations  are  valid  for  a  monitoring  surface  of  arbitrary 
shape.  Implementation  requires  the  computation  of  Green's 
functions  appropriate  for  the  application  of  interest. 

A. 3  FORMULATION  FOR  SURFACE  WAVES 

For  many  problems  the  primary  interest  is  in  the 
normal  mode  contribution  to  the  solution  for  waves  in 
layered  media.  For  such  problems  the  analysis  is  more 
conveniently  done  in  the  frequency  -  wave  number  domain. 

We  will  develop  the  solution  using  Green's  functions 
computed  following  Harkrider  (1964,  1970),  though  the 
theory  is  easily  modified  to  use  alternative  techniques 
for  computing  Rayleigh  wave  modes. 

For  the  azimuthally  averaged  Rayleigh  wave  displace¬ 
ment,  Green's  functions  for  a  particular  mode  are: 


i7T  ~R  eo  X(Z)  Jl(kC) 


i7r  £r  o  W(z)  Jo(Kr) 
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h{2) (kr) , 


(A. 11) 
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where  X(z)  and  W(z)  are  the  normalized  horizontal  and  ver¬ 
tical  displacement  eigenfunctions  and  eo  is  the  ellip- 
ticity.  The  AR  is  the  amplitude  response,  k  is  wave 
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number,  r  is  the  radial  distance  to  a  point  on  the  moni¬ 
toring  surface  and  R  is  the  range. 

The  stress  Green's  functions  are  easily  computed 
from  (A. 10)  using  the  following  identities  given  by  Hark 
rider  (1970): 


(A. 12) 

3W(z)  _  .  ( AX(z)  +  U»).\ 

3Z  *  \  A  +  2y  j 

where  t(z)  and  £(z)  are  the  normalized  horizontal  and  ver¬ 
tical  stress  eigenfunctions. 

Anelastic  attenuation  is  incorporated  using  the 
standara  e"Y(w)r  operator.  The  parameter  Y  is  computed 
from  specified  values  of  Q  for  each  layer  in  the  crustal 
model . 
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